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Abstract 



A simple model is considered to study the effects of finite size and internal 
structure in the tunnelling of bound two-body systems through a potential 
barrier. It is demonstrated that these effects are able to increase the tunnelling 
probability. Applications may include nuclear fusion, hydrogen atom and 
Cooper pair tunnelling. 
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I. INTRODUCTION 



The processes using the mechanism of quantum tunnelhng are abundant in different 
areas of physics. Chemical reactions at low temperatures, diverse phenomena in sohd state 
physics, nuclear alpha and cluster decays, fission and fusion processes in thermonuclear 
reactors, shortly after the Big Bang and inside stellar matter proceed through the Coulomb 
barrier. The conventional textbook treatment is, as a rule, limited by considering tunnelling 
of a point-like object. It is known that the intrinsic degrees of freedom of a system modelled 
by a barrier are important, and one needs to take into account the probabilistic distribution 
of barriers [1] in order to explain the subbarrier fusion. The friction-like processes under 
the barrier hinder the tunnelling [2]. At the same time, interaction of a charged particle 
with virtual electromagnetic field (an analog of the Lamb shift in bound states) increases 
the tunnelling probability [3]. In reality, in many cases one has to deal with the tunnelling 
of a composite object with its own intrinsic degrees of freedom. 

An important question of how the tunnelling reactions are influenced by the finite size 
and the structure of the tunnelling object remains poorly understood. For example, the 
dcuteron is a weekly bound complex with a size that is large on nuclear scale. How does 
this size influence the fusion reaction d + d — >^He deeply under the Coulomb barrier? There 
are no experimental laboratory data for the fusion cross section at low (~ keV) deuteron 
energies, since the cross section is too small. However, there is a discrepancy between 
the standard theory and the existing measurements for rfrf-rcactions in solids at energy 
between 10 and 2 keV [4,5]. This is an area of extreme importance for astrophysical and 
thermonuclear applications, especially with the perspective of radioactive beam facilities of 
new generation. Currently, the nuclear cross sections for low energies are determined by 
extrapolation based on the results for point-like tunnelling particles. 

Another example of great interest in relation to various tunnelling devices is given by 
electron Cooper pairs in superconductors. In "usual" superconductors, the size of the pair is 
so large that the electron binding is of minor significance for such processes, and the Joseph- 
son tunnelling matrix element can be taken as a product of penetrabilities for independent 
electrons. In high-temperature superconductors, the size of the pair is much smaller, and 
the binding may produce corrections. The tunneling of a bound electron pair through the 
potential barrier created at the "point contact" in semiconductors is also of interest because 
this process may be related to the mysterious 0.7(2e^//?.) structure in the measurements of 
conductance quantization in one-dimensional systems [6] (note that a bound pair appears in 
a one-dimensional system even for an arbitrary weak attraction) . One can also be interested 
in finite size effects for hydrogen tunnelling (a system of one heavy particle and one light 
particle). 

To provide an analytical estimate of the finite size effect we use a simple model. A 
system consists of two particles bound by the oscillator potential. An additional parabolic 
potential of positive or negative sign acts on one of the particles. This model allows one to 
study the limiting cases and shows that the finite size of a tunnelling system may be very 
important. Similar results for a pair of particles bound in an infinite potential box, while one 
of the particles encounters a rectangular barrier, were found in the old paper by Zakhariev 
and Sokolov [7], with the conclusion that the effective barrier is modified which facilitates 
the penetration. We show also that the electrostatic polarizabihty of a composite particle 
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adiabatically approaching the Coulomb barrier increases the tunneUing probabihty. 



II. TWO OSCILLATORS 



A. Normal modes 



We consider two particles interacting via the harmonic potential and moving in the 
external field that acts on one particle only. The Hamiltonian of such a system is 



H 



i + i + lk{x - yf + -qx\ 



(1) 



Here coordinates, x and y, and corresponding momenta, and Py, refer to a "proton" and 
"neutron", respectively; the potential (l/2)A;(a; — y)"^ describes the binding of the particles, 
and the potential {l/2)qx'^ acts on the "proton" only. The case oi q < models a Coulomb 
potential barrier for the "proton" in the process of the fusion of two "deuterons" into ^He 
or fusion of the "deuteron" with any other nucleus. 
The secular equation for normal modes reads 



4 2 
LU — U! 



'k 



kq 



+ — + 



0, 



(2) 



where the reduced mass = nixmy/ {rrix + friy). The solution may be conveniently written 
in terms of the parameters of bare frequencies and the ratio of the masses. 



A 



TTlx + m. 



(3) 



y 



Then the frequencies of two normal modes are 

o 1 



a;. 



at < 0, a;_ is imaginary. The case of equal masses, rUx — my = m, isA = l/2, when 
. 1 



1/2 + Q2 ± = —\q + 2k± Jq^ + 4A;2 

J 2m L * 



If the x-particle ("proton") is very light, mx ^ 0, A ^ 0, 



2 2 I o2 2 



0; 



(4) 



(5) 



(6) 



if the y-particle ("neutron") is very light, mx — > 0, A — > 1, two bare frequencies are normal 
modes. 



2 2 2 r>2 

UJ^ — V , UJ_ — \l. 



(7) 
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B. Wave function 



To find the spatial structure of tlie normal modes, we carry out the standard diagonal- 



ization of the Hamiltonian. First we make the scale transformation to new coordinates u 
and new momenta 



Then the Hamiltonian takes the form 



k f Un 



/m. 
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(8) 



(9) 



Since the kinetic form is invariant, we diagonalize the potential form by the orthogonal 
transformation 



Ux — Ci cos + ^2 sin (f), Uy — —^1 sin + ^2 cos 0. 



The transformed potential energy is 



where 



U22^k 



C0S( 



+ 



sm I 



sin cos ( 



H cos 0, 



H sm ( 



, I COS0 sin0 \ / sin0 cosi 
C/i2 = K I ^= + ^= I I ^= — 



/m 



y , 



/m,r 



cos sin (f). 



The diagonalization condition t/12 = determines the mixing angle 

2k^m^my 2u'^^X{l - A) 



tan(20) = 



Prom here 



qruy + k{my — nix) + ^^^(1 — 2A) 

2 n 



sin (20) 
cos(20) = 



2;y^y'A(l-A) 



^Q4 + i/4 + 2i/2Q2(i_2A)' 

02 + ^y2(i _ 2A) 



+ 1/4 + 21/2^2(1 _2A)' 
Therefore we get, as it should be, in agreement with (4), 

The wave function of the system is the product of two oscillator functions, 

^{x,y) = V^l(Cl;a;+)^2(6;^-)• 



(10) 
(11) 
(12) 
(13) 
(14) 

(15) 

(16) 
(17) 

(18) 
(19) 
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III. PHYSICAL EFFECTS 



A. Probability of tunnelling 

The most interesting case corresponds to negative q and, whence, negative a;^. The 
transmission coefficient for the parabolic barrier is given by [8] 

T=^, (20) 

where 

is the scmiclassical penetrability for a "one-way" traverse. Here E is the energy of the 
"deuteron" for motion along the coordinate ,^2, corresponding to the imaginary frequency; 
E = Q corresponds to the top of the barrier. The tunnelling probability is exponentially 
small for £■ < 0. A realistic case corresponds to A; a "deuteron" bound much stronger 

than Coulomb energy. In this case (for equal masses) 









V 2m I 



The factor 1 + \q\/{^k) provides an exponential enhancement of the tunnelling probability 
for the finite size "deuteron" in comparison with a particle of a zero size {k — 00) , 



where 

For the real deuteron this effect is too small but it increases if the "charged" particle is very 
hght, m^/niy -C 1, 



l + i^l- (25) 



^ vHx + my \ 8k rux 



B. Polar izability 

There is another effect that may be even more important than the modification of the 
barrier (change of cc;_) considered in the preceding subsection and analogous to the effect of 
ref. [7]. This is energy transfer from internal motion of the "deuteron" to the translational 
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motion which happens for a non-parabohc barrier. Let us consider a more reahstic model, 
where the potential barrier is smoothly goes to zero at infinity being parabolic near the top. 
The change of internal motion of the "deuteron" in the ground state occurs adiabatically 
as a transition from the unperturbed value hv/2 far away from the barrier to This 
energy difference is converted into energy of translational motion of the "deuteron" : 



3/2 



E = E(c.) + --^^E(^) + ^^y-^^ . (26) 

This correction also grows if the "charged" particle is light. 

For the real deuteron (or any charged system) this correction may be presented as a 
result of the deuteron polarization, 

E = E{oo) + -aE"^, (27) 

where a is the deuteron polarizability and £ the electric field produced by the Coulomb 
field. The internal ground state energy decreases by this amount, 5E = —aS'^/2, that 
is transferred into translational motion. Indeed, the electric field polarizes the deuteron 
(the electric dipole moment is d — aS). This creates an additional force, {d/dx)aS'^/2. 
Integration of this force gives a positive correction to kinetic energy, or, to present it more 
conveniently, a negative correction to the Coulomb barrier potential, 

1 Z2p2 

SU = --aS^ = (28) 

2 2r4 ^ ^ 

Again, this correction increases the barrier penetration factor. For the deuteron, the 
empirical value [9] is q;=0.70(5) fm^. We can estimate the change of the Coulomb barrier 
penetration factor due to this potential {U — Ze^/r + SU). In the WKB semiclassical 
approximation, 

where v is the relative velocity, m the reduced mass for two nuclei, and the cut-off radius 
where the validity of the potential approximated hy U = Ze'^/r + 6U is violated (rc ~ 
where is the deuteron radius). For the fusion of two deuterons, the correction due to 
the polarization potential is small (~ 1%). However, it rapidly increases with the nuclear 
charge. 



C. "Cold fusion" and muon catalysis 

In a process of muon catalysis, the fusion proceeds from a ground state of a muonic 
molecule, dd/i or dt/i. This process occurs due to zero-point oscillations in the ground state 
of the molecule, as soon as a neutron from one nucleus reaches the strong potential well 
of another nucleus. In the two-oscillator model with q > 0, the expectation values of the 
normal coordinates in the ground states are 
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Going back to the original coordinates, eqs. (8) and (10), we find 

h / cos^ (j) ^ sin^ 
2mx \ UJ+ uj- 



(x^) = -^f^ + ^V (31) 



,9, h f cos^ 6 sin^ (b\ 

^ = ^ ^ + ' 32 

2my \ ou^ ^+ J 

{xy) = 'r^^ ""'^ f- - -V (33) 

The osciUations of the "neutron" are described by {y"^). In the most interesting case k ^ q 
we have 

2\ ^ , R 



In this approximation the result for the "proton" oscillations (x^) is the same. Again, the 
finite size of the "deuteron" leads to stretching of the zero-point oscillations by a factor 



(1 + ^Jq/Ak) and to the enhancement of the fusion probability. 

D. Tunnelling of a bound electron (Cooper) pair 

The Hamiltonian of such a system may be modelled as 

H = + hix - yf + \q{x' + y'). (35) 

Then the frequencies of two normal modes are 

2 2k + q 2 (1 (oa\ 
m m 

At g < 0, a;_ is imaginary. It does not depend on /c, therefore the parabolic barrier for 
center-of-mass motion is not modified by the finite size of the system. However, there is the 
energy transfer from the internal motion to the translational motion if the potential barrier 
smoothly goes to zero at large distances being parabolic near the top: 



hJ2k/m nJ{2k + q)lm 
E = E{^) + ^ ^ (37) 

This energy transfer increases the tunnelling probability in the case of a smooth barrier, e.g. 
at the the point contact. The case of the tunneUing through the Josephson barrier is more 
complicated since the adiabatic approximation is not valid there. The situation might be 
closer to the regime of a sharp perturbation that deserves a special consideration. 
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IV. CONCLUSION 



Our simple model helps to understand that, even in a two-particle system, the tunneUing 
is a complex process that can be noticeably influenced by the flnite size of a tunnelling object 
and by its intrinsic degrees of freedom. We pointed out the necessity to study the behavior of 
the normal modes of a tunnelling system and its adiabatic polarizability (or deformability) 
by the barrier potential. The role of restructuring of complex particles in the process of a 
nuclear reaction was emphasized in a different context in refs. [10,11]. 
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the discussions of the dcutcron polarizability with W. Lynch. The work was supported by 
the NSF, grants PHY-0070911 and PHY-0244453, and by the Austrahan Research Council. 
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